Competing order, Fermi surface reconstruction, and quantum oscillations in 
underdoped high temperature superconductors 
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We consider incommensurate order parameters for electrons on a square lattice which reduce 
to d-density wave order when the ordering wavevector Q is close to Qo — {n/a,n/a), a being the 
lattice spacing, and describe the associated charge and current distributions within a single-harmonic 
approximation that conserves current to lowest order. Such incommensurate orders can arise at the 
mean-field level in extended Hubbard models, but the main goal here is to explore thoroughly the 
consequences within a Hartree-Fock approximation. We find that Fermi surface reconstruction in 
the underdoped regime can correctly capture the phenomenology of the recent quantum oscillation 
experiments that suggest incommensurate order, in particular the de Haas-van Alphen oscillations 
of the magnetization in high fields and very low temperatures in presumably the mixed state of these 
superconductors. For 10% hole doping in YBa2Cu306+«, we find in addition to the main frequency 
around 530 T arising from the electron pocket and a hole frequency at around 1650 T, a new low 
frequency from a smaller hole pocket at 250 T for which there are some indications that require 
further investigations. The oscillation corresponding to the electron pocket will be further split due 
to bilayer coupling but the splitting is sufficiently small to require more refined measurements. The 
truly incommensurate d-density wave breaks both time reversal and inversion but the product of 
these two symmetry operations is preserved. The resulting Fermi surface splits into spin up and 
spin down sectors that are inversion conjugates. Each of the spin sectors results in a band structure 
that violates refiection symmetry, which can be determined in spin and angle resolved photoemission 
spectroscopy. For those experiments such as the current photoemission experiments or the quantum 
oscillation measurements that cannot resolve the spin components the bands will appear to be 
symmetric because of the equal mixture of the two spin sectors. There is some similarity of our 
results with the spiral spin density wave order which, as pointed out by Overhauser, also breaks time 
reversal and inversion. Calculations corresponding to higher order commensuration produces results 
similar to anti-phase spin stripes, but appear to us to be an unlikely explanation of the experiments. 
The analysis of the Gorkov equation in the mixed state shows that the oscillation frequencies are 
unshifted from the putative normal state and the additional Dingle factor arising from the presence 
of the mixed state can provide a subtle distinction between the spiral spin density wave and the 
d-density wave, although this is very difficult to establish precisely. 
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I. INTRODUCTION 



Despite more than twenty years of intense effort, the 
telltale evidence of competing order parameters in the 
phase diagram of high temperature cuprate superconduc- 
tors remains an enigma. Other than antiferromagnetism 
of the undoped materials,'!' and stripe order at special 
doping and special materials,'^] not much is truly under- 
stood. Yet, the attractive option of classifying the pseu- 
dogap as a bro ken symmetry state has spurred intense 
interest .'SEEEEl Indeed, this option has the great power 
to unify the diverse phenomenology of these materials. 
We are strongly motiv ated by th e recent quantum oscil- 
lation measurements^ 9|io|ii|i2 [i3] ^^ address this question 
anew. These are the Shubnikov-de Haas (SdH) effect, the 
de Haas-van Alphen (dHvA)effect and the oscillations of 
the Hall coefficient (Rh)- The rapidly evolving experi- 
mental situation suggests, under various circumstances, 
evidence of commensurate and incommensurate order. It 
is perhaps best at this time to explore as thoroughly as 
possible both kinds of order within a very general frame- 
work. Such a framework is provided by the Hartree-Fock 



theory. 

Our basic assumption, justified when compared to ex- 
periments, is that the observed quantum oscillations in 
high magnetic fields and low temperatures in underdoped 
YBa2Cu306+5 arise from the Landau levels of the quasi- 
particles of the Fermi liquid corresponding to the puta- 
tive normal state. It is known, and shown here, that in 
the mixed state of many superconductors the oscillations 
remain unshifted in their frequencies from the normal 
state but exhibit increased damping. Thus, the beauty 
of the present experiments is that one can in principle de- 
termine the closed Fermi surfaces, through the Onsager 
relation that involves only the fundamental constants and 
the Fermi surface area, of the normal state without hav- 
ing to cross the upper critical field. At face value these 
observations have revealed some striking facts not no- 
ticed previously.'^ The most important of which is the 
existence of both positively and negatively charged quasi- 
particles in a hole doped cuprate at variance with the 
conventional wisdom. That these particles behave like 
quasiparticlcs in a Fermi liquid, even in the underdoped 
regime, which has so far been plagued with complex the- 
oretical concerns involving the proximity to a Mott in- 



sulatorflS is striking to say the least and should be of 
great value in uncovering the mystery of these enigmatic 
superconductors. One of us has elaborated on the strik- 
ing nature of these discoveries and has explained that 
the root of these observations is the Fermi surface re- 
construction due to broken symmetries P^ One of us has 
also elucidated the role of commensurate singlet d-density 
wave (DP W) to explain the oscillations of the Hall coeffi- 
cient .1^ The Hall measurement is more telling than either 
SdH or dHvA because it contains the striking informa- 
tion regarding the sign of the charge carriers. There is 
evidence of the presence of electron pockets in what is 
really a hole doped material, which is difficult to explain 
without spontaneously breaking symmetries in the puta- 
tive normal state, in particular by the featureless liquid 
of resonating valence bond theories.l^l 

For a particle-hole condensate, which DDW is, the 
symmetry of the orbital wave function does not constrain 
the symmetry of the spin wave function. ^^ Triplet DDW 
has not yet been adequately explored and remains a sub- 
ject for future work. The more recent experiments have, 
however, revealed some indication of an incommensurate 
order ing ,ISI which was only briefly touched upon previ- 
Q^g]^y|i6ii9| jjej-g -^g focus on the incommensurate case. 

Before we begin, it is important to state what a 
Hartree-Fock theory can or cannot accomplish. A 
Hartree-Fock theory can identify the possible broken 
symmetries and phases, but cannot establish that such 
symmetries are indeed broken. The prime deficiency is 
the absence of all fluctuation effects, but often such defl- 
ciencies can be remedied by educated guesses. Of course, 
deep inside a broken symmetry phase critical fluctua- 
tions are absent and some of the collective modes can 
be physically identified from the symmetries of the or- 
der parameters. Most importantly, symmetries protect 
the Hartree-Fock calculations; those properties that are 
determined by the symmetries alone can be determined 
in the weak interaction limit, where the calculations are 
better controlled, from a suitable effective Hamiltonian. 
The results protected by symmetries should be at least 
qualitatively valid even in the strong interaction limit .1^ 
However, the thermal properties in low dimensional sys- 
tems are badly predicted by the Hartree-Fock theories, 
because fluctuations are often very important. Nonethe- 
less, it is expected that the results at zero temperature 
retain a considerable degree of validity, except close to 
quantum critical points where quantum fluctuations be- 
come important. 

When a commensurate one-dimensional (1-D) charge 
density wave (CDW) is doped, the resulting charges may 
be viewed as defects in the CDW. If the interactions 
between them are sufficiently strong compared to their 
effective kinetic energy, then they will form an ordered 
lattice (which can be stabilized by a crossover to three di- 
mensions if there are many ID chains coupled together) 
of defects and the CDW will become incommensurate 
with the lattice. This can be energetically more favor- 
able than simply doping holes into the rigid band formed 



in the presence of CDW order at a fixed wavevector. The 
reason, in the latter case, is the single-particle energy gap 
of the CDW. In a two-dimensional DDW state, however, 
there are nodes in the order parameter, so it may be more 
favorable for doped holes to simply go to the nodes and 
expand into Fermi pockets (in quasi-lZ) systems, how- 
ever, this might not be possible, so doping must lead to 
incommensurability^^). However, this cannot persist in- 
definitely if the DDW state is stabilized by approximate 
nesting, since the Fermi surface would eventually move 
away from the nesting wavevector. Thus, incommensu- 
rate DDW order is a possibility, at least over a range 
of doping. In this paper, we explore the phenomenology 
of this possibility. But it is very difficult to determine 
if this incommensuration happens before the commensu- 
rate gap collapses and some other order takes over. In 
this respect the tendency of the spin density wave (SDW) 
to incommensurate is stronger because it is gapped every- 
where on the Fermi surface, being intrinsically a s-wave 
object. 

Near half-filling, the Hubbard model, the t — J model, 
and generalizations of these appear to have many phases 
which are close in energy. Thus, small perturbations can 
strongly influence the competition between them. Con- 
sequently, it should not be too surprising that experi- 
ments on the cuprate superconductors have also uncov- 
ered some evidence for a cornucopia of phases, particu- 
larly on the underdoped side of the phase diagram, which 
appear in particular materials and for certain values of 
the doping level, temperature, magnetic field, etc. De- 
pending on the interfacial energies between these phases, 
one way in which their competition can be resolved is 
through the formation of stripes^^ or other inhomoge- 
neous patterns of microscale phase-separation. Further- 
more, as the cuprates are doped, their Fermi surfaces 
evolve away from the two-fold commensurate nesting 
wavevector Qo = (7r/a,7r/a), so we would expect trans- 
lational symmetry-breaking order parameters to occur at 
incommensurate wavevectors. Hence, it would seem im- 
portant for any theory of the pseudogap to incorporate 
the tendency towards incommenuration, which we take 
as further impetus to study incommensurate order pa- 
rameters related to DDW order, the IC-DDW. Direct 
attempts to measure DDW ord er thr ough neutron scat- 
tering have neither ruled it out ESMl j^or unambiguously 
verified its presence, b ut th ere are intriguing suggestions 
that it may be presenlpSMl. 

In Sec. II the construction of both commensurate and 
incommensurate order parameters are discussed. This 
discussion serves to motivate our principal approxima- 
tions. Sec. Ill on fermiology is the heart of our paper. 
We begin with the calculation of the electronic spectra 
of a chosen IC-DDW order and deduce the correspond- 
ing areas of Fermi pockets that, through the Onsager 
relation, lead to the oscillation frequencies. The simi- 
larity of these results with the spiral SDW scenario is 
explained. Higher order commensurate DDW results are 
also discussed. These also turn out to be similar to the 



anti-phase spin stripe calculations. Thus, simply from a 
Hartree-Fock calculation it is not possible to single out a 
mechanism. One must invoke other considerations to dis- 
tinguish between mechanisms. Next we show that inter- 
layer coupling can lead to bilayer splitting of the areas of 
the electron pockets leaving the hole pockets essentially 
untouched. The subtleties of this result are discussed in 
some depth. Sec. IV contains a calculation of the self 
energy of the quasiparticles in the mixed state patterned 
after the work of Stephen,-2S- which we find to be by far 
the soundest approach, despite its deficiencies. The cal- 
culation for the nodal quasiparticles of DDW is a little 
subtle and leads to a vortex contribution of the scattering 
rate that we find very reasonable. If this calculation can 
be better controlled, it may be possible to distinguish 
between the SDW and the DDW scenarios. Given the 
approximate nature of this calculation, which can only 
yield a semiquantitative estimate, it seemed sufficient to 
illustrate our calculation with the commensurate case. 
Of course, the vortex contribution to the scattering rate 
is only a part of the total that also contains an impu- 
rity contribution. In the concluding section. Sec. V, we 
discuss our thoughts on microscopic models and whether 
or not incommensuration does actually occur. We also 
elaborate on experimental consequences of our work and 
point out the unresolved issues and future directions of 
research. 



the spin wave function. Thus, there can also be a triplet 
DDW which corresponds to circulating staggered spin 
currents,!^ as opposed to charge currents. If we set W{]i) 
by a function that transforms as the identity and remove 
the factor i, the order parameter would correspond to 
the conventional CDW. Similarly, the triplet DDW is the 
£ = 2 generalization of the conventional spin density wave 
(SDW) and is given by 
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a,k) = *W^k n • (T„,a'<5k',k+Qo , 



(3) 



where cr's represent the Pauli matrices, and n is a unit 
vector in the spin space. Once again, if we replace iW^ 
by unity, this order parameter is the conventional SDW. 
Because Wq is real, the order in Eq. ^ breaks time- 
reversal and hence represents a state with non-zero or- 
bital current configuration. This is most easily seen in 
real space, where Eq. (IT]) can be written as (the spin 
indices a suppressed for clarity): 



(4cr) = ±z(-l)^ 



hn'W^O 



Vr' 



R '.R 



(4) 



The d-wave hopping amplitude factor Vr'_r is 



Vr',R — —S^\R+ia — (^R'.R-xa + (^R'.R+ya + (^R'.R-ya 

(5) 
and also represents a current-conserving vertex centered 
at R as in Figll] (a) . A positive/negative sign in front of 



II. ORDER PARAMETERS 



A. Commensurate 

Before addressing the IC-DDW order parameter, let 
us STimmarize some of the features of the commensurate 
casc.^^ ""^ A commensurate singlet particle-hole conden- 
sate is defined by the order parameter 
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'C„ 



±iWk5a',Q<5k',k+Qo, 



(1) 



where Qo = (7r/a,7r/a). Although the right hand side 
should correctly contain a factor ±iWk, we will choose 
iWk for clarity, as no confusion should arise. The elec- 
tron destruction operator, Cka, is indexed by the wave 
vector k and spin a; a is the lattice spacing. We say that 
the order is commensurate, or strictly 2-fold commen- 
surate, because 2Qo is a reciprocal lattice vector. The 
form factor Wk transforms non-trivially under the point 
group of the two dimensional square lattice, defined by 
R = max + nay, where (to, n) is a set of integers positive 
or negative. It is related to the orbital angular momen- 
tum of the condensate. For DDW, 



Ww 



Wo 



(cos k^a — cos kytt), 



(2) 



which corresponds to angular momentum i — 2 and 
dxi_yi symmetry. For a particle-hole condensate, the 
symmetry of the orbital wave function does not constrain 




FIG. 1: The 6-vertex model generalization of DDW. The ver- 
tices (c)-(f ) are the singlet p-wave density wave vertices. Such 
a model has a thermodynamic phase transition to the low 
temperature DDW state but with no specific heat anomaly. 

each term in Eq. (|5| corresponds to current going out/in 
of the vertex. This analogy motivated a 6-vertex descrip- 
tion of the transition to a commensurate DDW state ,'^21 
where in addition to the d-wave vertices in Fig. [ll(a)-(b) 
one considers the effects of the p-wave vertices in Fig. [l] 



(c)-(f): Such a description was shown to have a finite 
temperature phase transition without any specific heat 
anomaly. As it will become evident below, one would 
have to include additional vertices to describe the tran- 
sition to a bond-incommensurate DDW state. 



B. Incommensurate 

To obtain an incommensurate version of DDW order, 
one first notes that since it breaks the Z2 time-reversal 
symmetry, one can construct antiphase-domain walls of 
mutually reversed current configurations. In Fig. [2] (a) 
we have an example of a bond-oriented domain wall. It 
can be written as. 



(4cr) = i{^iy 

where 



-"' ^Vr,,r X Q{ma) + zV^^,%, (6) 
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FIG. 2; (a) A bond-oriented domain wall, (b) The additional 
vertices ensuring current conservation at the domain bound- 
ary. 



is an overall antiphase modulation. The additional ver- 
tices V^,^^ in Eq. (pi are necessary to ensure current 
conservation along the domain wall boundary where the 
current is half of that in the bulk. Explicitly, they are 
given by the following: 



rCdgC 



Vr'3 — Vr/R^to^O + Vr, R(5„i,l 

Vr'^R = <5R',R-xa — ;^'JR',R+ya ^ ;TOR',R-ya 



(8) 



Vr'.R — ^R',R+Aa ~ ;:,^R',R+ya 



h 



R',R-ya 



and where V^, j^ and V§^, ^ are shown in Fig. (b). 

Note that the vertices in Eqs. (IS]) violate both d-wave 
and p-wave symmetries and hence a transition to bond 
IC-DDW order will not be in the universality class of the 
6-vertex model. On the other hand, a diagonal domain 
wall, as in Fig. |3] does not require extra vertices, and 
hence the transition to diagonal IC-DDW does belong to 
the 6-vertex model universality class. 

Because the contribution of V^/^^ is important only 
near a sharp bond-domain wall boundary, we drop it from 
now on. We will shortly see that at least in the single har- 
monic approximation, this is justified, for it leads to cur- 
rent conservation up to quadratic order in the inverse do- 
main wall spacing. Specifically, we can consider an array 
of antiphase-domain walls at positions ±d/2, ±3d/2, .... 
Because the array is periodic, we can focus on the region 
—d < ma < d, where due to the two domains at ±d/2 
the antiphase modulation 0(?7ia) in (pi is replaced by: 



6{ma+ -) — I 



with Ai = A/irl. 



d. __ v-^ , .nlma. ,„, 
^)-l= E Aicos{—-) (9) 
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FIG. 3: A diagonal domain wall. It can be entirely con- 
structed from the six vertices in Figjl] 
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Dropping the short range conributions Vj^, j^ and 
keeping the first harmonicP^lMI ; = 1 jn Eq. (lol amounts 
to considering an order parameter of the form: 



%'Cr; 



*(-I)"'+"'^cos(q-R)VR, 



R 



(10) 



where the incommensurability wavevector is proportional 
to the inverse domain wall seperation, so that for the 
bond-oriented domain in Eq. (|6l, q = (Tr/d, 0). Note that 
we regain the commensurate case of Eq. (H) as q — > 
since all domains are pushed to infinity in this limit. 



With the expression in Eq. (10) we can therefore ana- 
lyze a continuous transition to an IC state. It remains to 
check that the boundary contributions V^,^^ are unim- 
portant, by verifying that current is conserved at each 
site to quadratic order in q. 

The current along a bond R and R+ as, where s — x 
or s = y, arising from the IC-DDW state, Eq. (10 1 , is 
given by 



^R, 



Wo 



R4 



m-\-n ^ 



[cos q • (R -I- as) -I- cos(q • R)] . (11) 



Therefore, the total current flowing out of a site vanishes 
up to quadratic order in q. 

In the momentum space, the single-harmonic expres- 
sion, Eq. (lOl , looks particularly simple: 



(Ck'Ck) 



iWk [^k'.k+Q + 5k',k-Q] 



(12) 



As the ordering wavevector Q — (tt/c, tt/u) + q in 
Eq. ( 12 1 becomes commensurate, the above order reduces 
to that in Eq. (111. Moreover, because DDW is odd with 
respect to rotations by tt/2 and transpositions about the 
(tt, tt) direction one can check that in addition to domain 
walls, a checkerboard pattern results from the following 
momentum-space representation: 

(cLck) = ^Wk [4',k+Q + <5k',k-Q] - [Q ^ 0(Q)] (13) 

where the operation 0{Q) can be either a transposition, 
or a rotaion by 7r/2. Physically, the checkerboard pattern 
in Eq. ( 13 1 corresponds to simply superimposing two do- 
main walls rotated by 7r/2 with respect to each other. In 
deriving Eq. ( 12 1 we have considered bond-oriented do- 



mains, but the construction is easily generalizable for the 
diagonal domain walls in Fig. [Sj 



III. FERMIOLOGY 
A. Band structure 

We need to choose the band structure of the unrecon- 
structed Fermi surface. Although this is not precisely 
known, it is generally accepted that it should contain at 
least a nearest-neighbor, a next-nearest neighbor and a 
third-neighbor matrix element in a tight binding Hamil- 
tonian.ESlThus, a sufficiently general form is 

ek — — 2i(cos k^a + cos kyu) + At' cos k^a cos kyU 

- 2t"(cos 2kxa + cos 2kya) (14) 

In the past one of us has chosen^^ the parameters to be 
t = 0.3 cV, t' = 0.3i, and t" = t'/9.0, and to be consis- 
tent we shall adhere to these, although other reasonable 
choices will not change our main conclusions in the least; 



we could even set t" = or set it to a larger number, 
such as t" = 0.16tW^ This ambiguity is understandable 
because the precise results axe not known, either em- 
pirically or theoretically, especially in the underdoped 
regime. A more serious problem is to correctly recon- 
cile the picture of the Fermi surface emerging from the 
quantum oscillation observations with the angle resolved 
photoemission spectroscopy (ARPES).'^^ Finally, we will 
have to choose a chemical potential, n, to attain the car- 
rier concentration we desire, and this we will mostly take 
to be about 10% of holes, simply because many quantum 
oscillation experiments are available for this doping level. 
However, we will comment on what our theory yields as 
the the doping increases or decreases. 



B. Strict incommensurability 

In this subsection we address the situation when strict 
incommensurability holds. For example, Q = Qo + q = 
(tt/o, 7r/a) — 7r(2r7, 0)/a, where rj is not a rational number. 
It is of course a slight abuse of terminology to add the 
adjective 'strict' but it dispels any possible sources of 
confusion. In this case each and every k and k -f- Q 
correspond to a distinct point in the full Brillouin zone 
(BZ) of the underlying crystal lattice, and the fermion 
operators are similarly distinct. When necessary all sums 
will be carried out over the full BZ of the lattice with the 
conventional reciprocal lattice vectors. 

Consider now the mean field DDW Hamiltoniaan given 

by 



where 



Xw 



and 



n^J2 ^Ic^k^k 



Ak-M iGk ^ \ 

C.C. Ek+Q - M 

Ek ^ /^ ^*^k 

V C.C. ek+Q' - m/ 



Gk = (T^k-Wk+Q)/2 
Gk = (W/k-W^k+Q')/2 



(15) 



(16) 



(17) 
(18) 



The four component 



* 



(^kT ' '^k+QT ' '^ki ' Ck+QM ) • 



spmor la 1!^ 
It is useful to elaborate 
on this mean field flamiltonian. For a given wave 
vector Q = (vr/a, 7r/a) -|- q, the inversion symmetry is 
broken, while the pure imaginary character of the order 
parameter (and the compex Hermitian character of the 
Hamiltonian) implies broken time reversal symmetry. 
Note that the product of inversion and time reversal 
must be preserved as a physical requirement. Because 
the order parameter is a singlet in the spin space, this 
implies a very a simple matrix structure of the Hartree- 
Fock DDW Hamiltonian. The inversion conjugate of this 



order parameter, Q' = {n/a^n/a) — q, must belong to 
the down spin sector if Q belongs to the up spin sector. 
To further clarify the form of this Hamiltonian, note 
that the use of the full BZ allows us to appropriately 
shift the wave vectors to use identities like 

II4,aCk+Q,a(iGk)+h. C. (19) 

k 

= X! Ckc<Ck'Q(iWk) [(^k'.k+Q + f^k'.k-q] + h. c. 

k,k' 

for any given spin direction a. 

The Fermi surfaces obtained from the spectra of 
Eq. ( 16 1 are shown in Fig. [4] The pockets proximate 
to (±7r/a, ±7r/a)/2 are hole pockets while those close 
to (7r/a, 0) and its symmetry equivalents are electron 
pockets. To obtain this figure we have chosen: Wq — 
0.112 eV, ?7 = Om^ fj. = -0.256 eV, and nn « 10%. 
The corresponding dHvA frequencies determined by the 
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FIG. 4: (Color online) Reconstructed Fermi surfaces for in- 
commensurate DDW order: (a) the up spin case with wave 
vector Q = iji/a, n/a) + q; (b) the down spin case with 
Q' = (7r/a, 11 /a) — q; (c) the combined up and down spin 
spectra. 



Onsager relation 
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A{eF). 



(20) 



where ^(ef) is the cross sectional area of a closed Fermi 
surface, are Fi = 526 T (electron pocket), F2 = 1670 T 
(large hole pocket), and F3 = 250 T (small hole pocket). 



We have chosen the parameters to fit the data,'^' but the 
existence and frequency of the small hole pocket is a pre- 
diction for which there appears to be very preliminary 
indications requiring further investigations and confirma- 
tion.^ 

It is important to specify how we arrive at the total 
hole doping of 10%. Fermi surface is a topological invari- 
ant.'^^ Even when quasiparticles behave anomalously, as 
in one-dimensional electronic systems, this surface is still 
defined by the same topological invariant. A break up 
of this surface, termed reconstruction, requires a global 
deformation in the topological sense, most likely a macro- 
scopic broken symmetry. The key here is Luttinger's sum 
rule. The most general form of this sum rule states that 
the particle density is twice (for two spin directions) the 
volume of the wave vector space in d-dimensions divided 
by (27r)'' over which the real part of the single particle 
Green function at the Fermi energy is positive, which 
applies even to Mott insulators.'^ This is easily applied 
here. There is one big hole pocket for a given spin di- 
rection within the reconstructed HZ. Considering the two 
spin directions producing identical pockets but reflected 
with respect to each other, the big hole pockets corre- 
spond to a fraction of carriers, Xbh, to be 



a;bh 



,Ah(eF) 
' (27r/a)2 ■ 



(21) 



in terms of carriers per Cu. Here A^h^ep) is the area of 
a big hole pocket. Note that the normalization in Eq. 
(23 1 is with respect to the area of the unreconstructed 



BZ. Similarly, there is only one electron pocket in the 
reconstructed BZ and therefore the fraction of carriers, 
cCe, is, considering both spin partners. 



= 2 



Acjep) 
(2^/a)2 



(22) 



But there is also one small hole pocket per reconstructed 
BZ and therefore 



^sh 



^Ashjip) 
' (27r/a)2 



The hole doping, n^, is 



nh 



•^bh H~ -^sh '^c 



(23) 



(24) 



It is important to compare the DDW case with the 
spiral SDW case. The mean field potential, T^pirai, is 

^spiral = ^0 [S^ cos(Q -11) + Sy sin(Q • R)] . (25) 

where S^ and Sy are the conventi onal s pin half opera- 
tors. As pointed out by Overhauser,H5Eil this interaction 
Hamiltonian breaks time reversal and inversion. It im- 
mediately follows that the mean field Hamiltonian matrix 
for the spiral SDW is 



Yv 



^ek-M A 

A Ek+Q - /-i 

Ek-M A 

A ek+Q' 



(26) 



flit,, 


3-Qt L 

;: ^ 




3- r 


^ ; 



-3-2-10 1 2 3 

(a) 







flit,. 




3 


Q'i. 






2 
1 


/^ 


7 


^ 











n 


ak^ 






u 


-1 
_2 


^^ 


:^ 


^ 




— 3 




rn 






-3 -2 - 


-10 12 3 





(b) 







flit,. 






3 

2 

1 



-1 

-2 

-3 




J 


fl/l. 




r 


^ 

1 




-3 -2 


-1 1 


2 3 





(0 



FIG. 5: (Color online) Reconstructed Fermi surafces for in- 
commensurate spiral SDW order: (a) the up spin case with 
wave vector Q = (t^/cl, it /a) + q; (b) the down spin case with 
Q' = (tt/o, 7r/a) — q; (c) the combined up and down spin spec- 
tra. Note that even though up and down spins are mixed we 
continue to label the spectra with the spin directions, since 
Q and Q' defines the cases unambiguously. 



where the 



("-kT' 



,t 



'^k-l-Qi'^ki 



four 



-k-l-Q'T 



component 
). Here, 



$ 



spmor 

A, the spiral SDW 
gap parameter is of course real. But it is trivial to see 
that the eigenvalue structure is identical to the IC-DDW 
because both inversion and time reversal are broken. In 
contrast to IC-DDW, however, the up and down sectors 



are mixed instead, as it should be, because DDW is a 
singlet in spin space, while SDW is a triplet, bearing 
in mind that both condensates are in the particle-hole 
channel. We can redefine, however, the spinors, such as 



C^k+Q.a 



*^k,Q 

«eQ/3Ck+Q,/3 



(27) 



similarly for Q ^ Q'. The antisymmetric tensor Cq,^ is 
defined by en = e|| = £22 = exi = 0, £12 = ej; = 1, 
and £21 = £41 — —1. With the redefined spinors the 
Hamiltonian matrix is now 



Y^ = 



y^£k - /^ lA 

-lA £k+Q - 



V 






Ek- M 

iA 






-iA 



(28) 



which closely resembles Eq. ( 16 1 



Ck+Q' — ^y 

The transformation 



clearly leaves the anticommutation rules unchanged. 



The Fermi surfaces resulting from Eq. ( 26 ) are shown 
in Fig. [5] They are essentially indistinguishable from 
Fig. |4] This is despite the fact that the SDW order pa- 
rameter is a triplet and is gapped, separating the valence 
band and the conduction band, over the entire Fermi sur- 
face, in contrast to the DDW order parameter. This is 
unfortunate because from dHvA frequencies, there is no 
way of telling if they arise from an incommensurate spi- 
ral SDW or an incommensurate DDW; one must invoke 
other considerations. In constructing this Fermi surface, 
we have chosen the following parameters: A = 0.08 eV, 
1] = 0.08, Ai = -0.27 eV, and Uh ~ 10%. The correspond- 
ing dHvA frequencies are Fi — 533 T (electron pocket), 
F2 = 1667 T (large hole pocket), and F3 = 274 T (small 
hole pocket). We suspect that a more diligent optimiza- 
tion of these parameters could be performed, but it will 
not change the main observations. 



C. Higher order commensurability 

As an example, consider higher order incommensura- 
tionofQ = Qo-f(2x|,0) = f(|,l). 



With the 8-component spinor defined by xj^ — (c^t 



„t 



as 



J 



^k+Q,Q' ^k-|-2Q,Q' 



n 



/ ,Xka^k,aXka 

k,CK 



"k-fso)' ^^"^ Hamiltonian can be written 

(29) 



Now the up and down spin sector eigenvalues merely duplicate each other, and we can consider simply one of them: 

*Gk 



/ Ck-M 
c.c 








\tG 



k+7Q 



Ck+Q - 

C.C 









M 





iGk+Q 
ek+2Q - 

c.c 







fJ- 






iGk-|-2Q 
£k-|-3Q — A^ 

c.c 











^Gk-fSQ 
ek-F4Q — /^ 

c.c 











iGk-|-4Q 
ek+5Q 

c.c 




M 










*Gk+5Q 
ek+6Q — A 

c.c 



— iGk+7Q \ 







iGk.+6Q 
fk-f 7Q ~ ^J' / 



(30) 



In the real space, we have 
(4'CR> - ^ E<4'Ck) exp [-i{k' • R' - k • R) 



k'k 



= ±^ ^ Gkexp [ik ■ (R - R')] exp [-iQ ■ R'] 



'N 
iW< 



±^(-l)"'+"'(vk^R + ^C/R^R 



(31) 



where R' — {m'a,n'a), and Vr'.r and C/r'^r are of the 

form;'*^ 



V, 



R',R — 



1 + COS 2Trr] 



(^R',R+ax + <5R',R-ax) 



(<5R',R+a.y + (^R'.R-ay) 

sin 27rry sin 2m'TTri 



cos 2m'nri 

(OR',R+ax ~ (JR'.R-ax), 



Ur> 



R 



1 + COS 2TrT] 



(<5R',R+ax + (^R'^R-ax) 



((^R'.RH 



ay 



JR',R-ay 



y) 



sin 2m'Trri 



sin 27777 cos 2m'TTri 



(<^R',R+ax — <^R',R-ax)- 

(32) 



We compute the current pattern from 

^R',R = «[(cLcr) - (C^^CR')] 
= -M/o(-l)"'+"'Vk',R 



(33) 



As an example, the current pattern with r/ = 1/8 is drawn 
in Fig. |6] The relative amplitudes of currents on bonds 
are labeled by different arrows. Because Qy = tt/o, the 
currents on a vertical slice have same amplitudes but al- 
ternating directions. Qx = i^ I a — 2'Kri/a — Sn/Aa results 
in a modulation of a period of 8 lattice spacing along the 
horizontal direction. 

From J2keBZ ^k = 0, it follows that (c|^cr) = 0. This 
result is in disagreement with Ref. 142 [ where the sum was 
incorrectly performed over the RBZ. Namely there is no 
site charge modulation in the higher order commensurate 
d- density wave state. However, note that ?7r',r ^ 0, for 
R ^ R', and there are bond charge modulations. 

The Fermi surface corresponding to the spectra of 
Eq. (30 1 is shown in Fig. It] It is not essentially dif- 
ferent from the mean field theory of 1/8 magnetic an- 
tiphase stripe order. ^^ This higher order commensura- 
tion generically produces complicated Fermi surfaces, in- 
volving open orbits, hole pockets and electron pockets. 
However, it is difficult to satisfy simultaneously the con- 
straints of the Luttinger sum rule, the periodicity of the 
oscillations, and the negative sign of the Hall coefficient. 



It is not clear how this picture can be consistent with 
experiments, unless disorder, fluctuations, or magnetic 
breakdown'*'^ conspire delicately to reproduce the exper- 
imental observations, which of course can not be ruled 
out. We shall not pursue this approach further. 
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FIG. 6: Current pattern for Q = (|^, ^). The relative magni- 
tudes of the currents are depicted by the arrows in the legend. 
Note the antiphase domain wall structure. 




FIG. 7: Reconstructed Fermi surafces for a plausible high 
order commensurate DDW with Q — ^(|,l)and Wo = 
0.0825 eV. There are electron pockets, hole pockets and open 
orbits. Note that the figure is depicted in the extended BZ 
for clarity. 



D. Interlayer tunneling and bilayer splitting 

Bilayer coupling, t±(k), in YBCO is parametrized in 
terms of a momentum conserving tunne ling matrix ele- 
ment. For tetragonal structure it i^nSEl 



i^(k) 



4 



[cos(fca;a) — cos(fcj,a)] 



(34) 



where a is the lattice spacing. The tunnehng Hamiltonian 
Hi2 is given by 



H 



12 



/^ 



i±(k)(c. 



t(l) (2) 
k,a k,a 



1 



2). 



(35) 



The superscripts on the ferniion operators refer to the 
layer index. For simplicity consider first the commen- 
surate case. The Hartree-Fock approximation, Hq — 
Hi + H2 can be written as the effective DDW Hamil- 
tonians 



kG-R-BZ,Q 



keRBZ,a 

+ (f ^2). 



(36) 



The reduced Brillouin zone (RBZ) is bounded by ky ± 
kx ~ ztw/a. With the choice of the quadratic Hamilto- 
nian Hq in Eq. (36), it can be easily diagonalized along 
with Hi2. This is a first order degenerate perturbation 
theory. Because t± will turn out to be so small, we do 
not expect a large correction. 

At each wave vector k in the RBZ, we need to diago- 
nalize a 4 x 4 matrix to extract the energy eigenvalues. 
This matrix ia^ 



iWu 







Ek iWk tj_(k) 

-iWw Ck+Qo t±{k+Qo) 

t±{k) Ek iWk 

t_L(k+Qo) -iWi, Ek+Qo 



(37) 



From t±{\i.) it is clear that the electron pockets will be 
much more affected by it than the hole pockets. The 
bilayer splitting of the electron pocket frequency in the 
dHvA meausurement, if it is to occur, should be smaller 
or of the order of the half-width at the half-maximum of 
the peak in the Fouri er s pectra, otherwise it would have 
been already resolvecpEl, Combined with the Luttinger 
sum rule this provides a strongl constraint on the chosen 
parameters. 

As an illustration for YBa2Cu306.5, we choose t± = 
8 meV, and Wq — 0.0825 eV, and the chemical potential 
fi is set to —0.2627 eV, which leads to a hole doping 
of Uh ~ 10%. The corresponding dHvA frequencies are 
Fi « 944 T, F2 « 967 T, F3 « 570 T, and F4 « 450 T. 
The frequencies Fi and F2, corresponding to the hole 
pockets are essentially unchanged to our accuracy, while 
F3 and F4 correspond to the electron pockets split by the 
bilayer coupling. The Fermi surfaces are shown in Fig. |8] 

The generalization to the incommensurate case is 
straightforward and has been discussed previously.'^ 
Here we do not pursue this further but turn to a number 
of important conceptual issues. 

1. The inter layer tunneling Hamiltonian is real and 
cannot lead to orbital currents flowing between the 
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FIG. 8: (Color online) Fermi surface split by bilayer cou- 
pling for commensurate DDW. The hole pockets centered at 
I (±7r/a, ±7r/a) are essentially unsplit by the bilayer coupling, 
In contrast, the electron pockets centered at (-Tr/a, 0) and sym- 
metry related points are split as described in the text. 



layers. This is actually quite fortunate because 
the DDW order parameters within the planes re- 
main undistorted. The sole effect of the interlayer 
Hamiltonian is to lead to linear superpositions of 
the independent quasiparticle states, which lead to 
a bilayer splitting, as in a generic two-state system. 

2. If we are to choose opposite phasing of the orbital 
currents of the two layers, and change the matrix 
to 



ek 


*VFk 


i±(k) 





-iW^ 


f=k+Qo 





ti(k+Qo) 


i±(k) 





Ek 


-iW^ 





i±(k+Qo) 


^W^ 


Ek+Qo 



(38) 

there is no bilayer splitting of the areas of the Fermi 
pockets. On the other hand, the ground state en- 
ergy, to a very high degree of precision, is identical 
to the previous cases even though the spectra are 
changed. We suspect that this equality is exact but 
have not been able to prove it analytically. There- 
fore, at the level of approximation discussed here, 
energetics do not allow us to choose between the 
alternatives and the correct phasing should be de- 
termined from experiments. 

We have argued that the interlayer tunneling 
Hamiltonian is both real and bilinear and therefore 
does not result in a flow of current between the 
layers. The only way such currents can appear is if 
there are further four fermion terms connecting the 
two layers such as particle-hole or particle-particle 
pair hopping terms. This would allow us to form 
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further condensates made out of the two layers and 
would seriously distort the planar DDW order pa- 
rameter. This has apparently gone unnoticed in 
Ref. HH The distortion of DDW order parameter 
must be thought through anew. 

There is a further important consideration that is 
worth spelling out carefully. One of us recently arguecpSl 
that the first order interalyer tunneling can be traded by 
a second order Hamiltonian in which virtual pair hopping 
of both particle-particle and particle- hole are present. 
The second order Hamiltonian involves four Fermi inter- 
actions from which appropriate condensates can be con- 
structed. It was shown there that for such a Hamiltonian 
antiferromagnetic arrangement of the orbital currents of 
the layers in a bilayer unit would be lower in energy than 
the ferromagnetic arrangement. There were two argu- 
ments behind this effective Hamiltonian. The first con- 
sisted in noticing that the DDW quasiparticle spectra 
are substantially gapped at the antinodal points. Alas, 
given the existence of electron pockets in the quantum 
oscillation experiments, this mechanism could not be op- 
erative because there is no longer any gap at the antin- 
odal points. The second argument involved non-Fermi 
liquid spectral function. One can easily see that with 
increased doping there would be Fermi surface reconnec- 
tions near optimal doping and the electron pockets will 
evaporate. But at such doping, close to the middle of the 
superconducting dome, there will be a quantum critical 
point and one can expect the non-Fermi liquid behav- 
ior.HZl The systematics of the superconducting transition 
temperature of optimally doped homologous series can 
still be correctly given by the second order pair hopping 
terms .'3Sj 



IV. QUANTUM OSCILLATIONS IN THE 
MIXED STATE 



There is considerable evidence that the quantum os- 
cillation frequencies remain unshifted from the normal 
state in the mixed state of a wide class of superconduc- 
tors, although there is increased dampingarising from 
vortices.l^SI We shall closely follow Stephen^^ where this 
problem is solved for a conventional metal in the nor- 
mal state and a s-wave superconductor with a gap Ag. 
Since the Bogoliubov quasiparticles do not couple mini- 
mally to the gauge field, they do not form Landau levels 
and clearly can not be the source of quantum oscilla- 
tions. ^° The quantum oscillations must then arise from 
normal quasiparticles. By solving the Gorkov equations 
in the mixed state Stephen arrived at a formula for the 
self energy, which for high Landau levels, relevant for the 
present problem, is 




where LOr 



Anrnhio, 
^ and m 



-msgniuj) ■ 



M 



hUJr 



(39) 



is the corresponding effective 
mass. The real part of the self-energy shifts the chemical 



potential and the Landau level positions equally. Hence, 
the real part of the self-energy does not affect the oscil- 
lation frequency. The imaginary part of the self-energy 
leads to the following scattering rate 



ttA^ 

V hhJr 



fi\hujc 



A^(l 



B 



B, 



c2 



\^JL\huJc 

(40) 
where the filling fraction is v — |//|/(/iWc) and Bc2 is the 
upper critical field. 

To illustrate our main points, we shall treat the simpler 
commensurate case. The extension to the incommensu- 
rate case is straightforward. There are two cases to con- 
sider. The effective mass corresponding to the electron 
pocket within the commensurate DDW theory, obtained 
by expanding around (tt/o, 0), is given by 



2t' + At" ^ I a^ 



2m 



(41) 



Note that t does not enter in the leading order. As- 
suming, Wo = 0.0825 eV, /x = -0.2642 eV, and Aq = 
10 mcV, we get 



1 



= 2.9 X 10^2 sec-\ 



(42) 



where we have used i? = 40 T and Bc2 = 60 T for the 
purpose of illustration. Using Eq. (41 1 and the parame- 
ters given above, we find that m* — 1.27fMei where m^ is 
the free electron mass. 

A more interesting situation arises for the hole pocket, 
which, to an excellent approximation, can be described 
by the nodal fermions of DDW. For the purpose of il- 
lustration, we choose again the commensurate case. The 
Stephen formula has to be rederived for nodal fermions. 

If the Fermi surface reconstruction is due to SDW, the 
excitations of both the electron and hole pockets can be 
approximated by the non-relativistic Schrodinger equa- 
tions with the appropriate effective masses and Stephen's 
formula will apply without further modifications. As an 
aside, we note that for the spiral SDW given in Ref.'F, the 
effective masses corresponding to both electron and hole 
pockets are considerably smaller (of the order of 0.5 rrie) 
and do not seem to agree with the experimental results. 
Similarly, the Dingle factors for the spiral SDW are differ- 
ent from those of the DDW order simply because SDW 
is gapped everywhere and the excitations are given by 
non-relativistic fermions. 



Commensurate DDW 

The linearized two component Hamiltonian for DDW 
quasiparticles, corresponding to the node (7r/2a,7r/2a), 
in the presence of an external magnetic field is given by 



/■ 



H= / dri}\r)H{v)il){r) 



(43) 
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where ^(r) is a two component spinor. The kernel H is 
given by 

H = vf{Px - eAx/c)a3 + voiPy - eAy/c)a2 (44) 

where cr's are the standard Pauh matrices; vp is the ve- 
locity orthogonal to the Fermi surface and vd is the veloc- 
ity tangential to it; A^ and Ay are the vector potentials 
and c is the velocity of light. The anisotropy between 
X and y directions can be removed with the redefinition 
of the coordinates x — y^vjj/vpx and y = ^/vp/voy- 
The measure of the integration remains unchanged, as 
dr = dr, and H becomes 

H = y/VFVD{Px - eAs:/c)a3 + y/VFVoiPy - eAy/c) (72 

(45) 
Under the unitary transformation ip —> U^, where U — 
(1 + iCT2)/\/2, H -^ H', such that 



H' 



= V^fvd{Px - eAi/c)ai + y/vpVD {Py ~ eAy/c)a2 

(46) 



Hence, the energy eigenvalues are 



e„,La = OL^/nfub 



ah\/2nvpvi) / 1 1 



(47) 



for n > I, where a = ±1 stands for the particle and 
hole branches of the spectra. Here a) — \/2vpVd/Ib, k — 
21:711/ Ly, m = 0,±1,±2,..., and the magnetic length 
Ib — yJhc/eB. The rescaled length Ly = ^Jvp/voLy. 
For n — 0, 



■^ofc 



eo = 0. 



(48) 



When sgn(ei3) > 0, the spinor wavefunctions in the Landau gauge A(f) = (0,55;, 0) are 






S(^S.fe(f),o) 



where 



'/5n,fe(r) = 



^lBLy2''n\ 



V2' 



exp —iky - 



2 ^2 



(x-kli,) 
21% 



H,, 



kf, 



and Hn are the Hermite polynomials. The matrix Green's function is defined as 

Go(ri,r2,«cj) = > — +> ^-—^ , 

■^ lU! — e„^Q -I- /i ^ ILU + ^ 

n>l.k,a k 



(49) 
(50) 

(51) 
(52) 



After performing the sum over k and a, the Green's function can be expressed in terms of redefined coordinates as 



Go(fi,f2,Ja;) 



2ttII 



~ir ^^P ( ^~ITr~ ) exp(i02i) 



4/2 



E 



iuj + fi 



Ln 



2/1 



n>l 






2n[{^u^ + m)2 - el] ^ S*2,Li_, (%f ) 





1 (I 
iu + ii\Q 



(53) 



where (j)2i = (^2 + a^i)(j'2 — yi)/2?^' -^21 = (2^2 — 2^1) + i{y2 — Vi), and L™ are the associate Laguerre polynomials. 
The Gorkov equation for the normal component of the Green's function of a s-wave superconductor is 

G(r,r',iw) = Go(r,r',itj) - / dridr2Go(r,ri,iw)As(ri)Go(r2,fi, -it^)A,*(f2)Go(f2,r', iw) (54) 

For a disordered configuration of vortices (for other choices, see StepherP^ 

(f2-fi)2 



T/(ri,f2) = (A,(fi)A:(f2)exp(2i02i)) = A.^exp 



2P 



(55) 
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From the Gorkov equation we can identify the real space self-energy for the average Green's function to be 

i;(ri,f2) = F(fi,f2)exp(2i02i)Go(f2,fi,-iw). 
The matrix elements of E(ri,r2) in the Landau level basis are given by 

^ni,fci,ai;«2,fc2,a.M = J ^^f ldr27/;^^ -^ ^^^ (f i)F(f i , f 2) exp(2z02l)Go(r2, ?! , -iw) V'„,^fe, ,„, (rj) • 

After performing the integrations we obtain 



^ni,ki,ai-n2,k2,a2^^^) ~ 9 ^ri^^^^'^kxM 



n=0 ^ ' " 






00 ^ 



n^O^^'^^ + ^^^^^^+l 



(56) 



(57) 



(58) 



r 



where 



{n + ni)\ 



When both n and ni are large and comparable to the 
filling fraction v ^ {\fi\/hCj)'^, 



1 



: exp 



4n 



(60) 



and E„(ia;) ^ E„j(iCLi). For large n and rii, using the 
asymptotic form of Im.n, and converting the the Landau 
level sums into integrals we find 



,(iw) 



Ay/n{hu}) 



:A2 



M 



(l-«i)(l 



- i sgn(cj/i) 
a2). (61) 



It should be noted that for only the hole branch of the 
spectra is the self-energy nonzero. When n ^ v , we can 
neglect the real part of the self-energy and find the vortex 
contribution to the scattering rate to be 



7rA| 

4 hij 



^VWk 



(62) 



where the effective mass for DDW quasiparticles and the 
cyclotron frequency analogous to the non-relativistic case 
are defined dimensionally as 



m 



VpVD 

eB 



(63) 
(64) 



Note that being relativistic Weyl fermions, there is no 
real mass associated with them. Here vp — 2\/2at/h; 
to leading order neither t' nor t" enter. Similarly, vjj — 
Woa/\/2h. Assuming, again, Aq = 10 meV, and Wq — 
0.0825 eV, we get 



1 



28 

Vb 



1 - 



B 



X lO^^sec-i 



(65) 



where B is in units of Tesla. For B = 
60 T, we get 

— = 1.5 X lO^^sec-^ 



40 T and Bc2 



(66) 



Despite the relativistic nodal fermionic character of the 
quasiparticles, the quantum oscillation formulas at zero 
temperature, using the method of Ando,^^ are formally 
identical to the non-relativistic case with appropriate re- 
definitions. The damped oscillatory factor for the funda- 
mental i^ 



-Tt/ll 



^" COS {2ttF/B) 



(67) 



where F is again given by the Onsager formula. If we 
also include the Dingle factor arising from the impurity 
scattering with a scattering time tj, we need to multiply 
the above formula by a factor e~'^/"<:'^'. At a temperature 
T, the amplitude. A, is 



Aoz 



sinh( 



2-R^kBT 
2Tr^kBT 

hut* 



(68) 



We emphasize that m* here is a parameter, not the 
effective mass of the nodal fermions, which are actu- 
ally massless. Thus, this formula could not be the 
conventional Lifshitz-Kosevich formula derived for non- 
relativistic fermions, despite its formal similarity. From 
Eq. ( 63 ) , the m* for the parameters given above is 
2.72 rrie, more than a factor of 2 larger than the non- 
relativistic mass corresponding to the carriers comprising 
the electron pocket. 

If we worked with a d-wave superconductor, then the 
Gorkov equation is more complicated due to the nonlo- 
cal gap function and there will be integration over four 
variables. However, the final answer for the scattering 
will not be significantly different from the s-wave answer. 
We can just replace A^ by the Fermi surface averaged 
(jAkP), which is the same as averaging over extremal 
orbit in 2D. 
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V. CONCLUSIONS 

A. Thoughts on microscopic models 

An effective Hamiltonian that can lead to both DDW 
and d-wave superconductivity (DSC) in Hartree-Fock 
theory is the followingP2l 

it — — / ^ tij [Cia^ja + h.C. I —tc 2_^ '^ia'^ja'^ja'^i'a' 

+ ;7^n,^n^-j^+l/^n,nj- (69) 

In this formula, ty is hopping matrix element with ty = t 
for nearest neighbors, t^ — —t' for next nearest neigh- 
bors, etc. The operator n^Q, is the density at a site i and 
spin a. On the other hand tc is a correlated hopping 
term which simultaneously hops an electron from site j 
to site i and hops an electron from i' into the vacated 
site j. The on-site and nearest-neighbor repulsions are, 
respectively, U and V. The repeated spin indices a, a' 
are assumed to be summed over. 

At the mean field level, the energetics of the DDW 
order is described by the effective interaction 



n 



DDW 



-5DDW / /(k)/(k')c[+Q„Cka4,„,Ck' + Qa', 
k,k' 

(70) 
with goDW = 24tc + 8V^ and /(k) = (cos fc^a — cos kyo). 
A similar reduced Hamiltonian for DSC has the coupling 
.9DSC — 12ic — ^V .^^ Note that the processes responsible 
for both DSC and DDW are essentially kinetic. 



A Hamiltonian given in Eq. ( 69 1 has a d-wave super- 



conducting ground state over a range of dopings and an 
antiferromagnetic ground state at half-filling. Since this 
is the sine qua non for any description of the cuprates, 
we believe that this is a good starting point for further 
calculations. Correlated hopping terms naturally arise 
even from the one-band Hubbard model away from half- 
filling, where their contribution can be computed in t/U 
perturbation theory: 



tc 



nh 



C/2 



(71) 



At dopings nh — 0.1, for t/U ~ 0.2, we would have 
^DDW = 2gDSC ~ 0-li ~ J/'^i J being the antiferro- 
magnetic exchange constant. Thus, not far away from 
half-filling, correlated hopping terms are certainly not 
negUgible. 

Preliminary calculations^* at the Hartree-Fock level 
shows that it is difficult to make DDW incommensurate 
in the relevant underdoped regime. The physical reason 
is that DDW has nodes. Our mean field results suggest 
that it costs less energy to dope holes by opening pockets 
around the nodes rather than by making DDW incom- 
mensurate, as long as there is sufficient density of states 



left at the antinodes to support commensurate DDW. 
To confirm the above intuition, we tested the effect of 
gapping the nodes by either introducing s-wave order or 
including a id^^y piece^to DDW. We could see that at fi- 
nite chemical potential, the effect of increasing the nodal 
gap created by the extra order is to induce a first order 
transition to an incommensurate state when the nodal 
gap amplitude is larger than a critical value, of the order 
of the dj;2_y2 amplitude itself. Of course Hartree-Fock 
approximation need not be the full answer because the 
notion of topological doping^ in the context of stripe 
order can be equally operative in the present context, a 
point that requires careful further work. 

It is amusing to note that spiral SDW is also generically 
energetically unfavorable compared to coUinear SDW in 
the weak coupling Hartree-Fock approximation, as was 
pointed out by Overhauser.*'^ The strong coupling anal- 
ysis can in principle lead to spirals, but these calculations 
are difficult to justify.!^ 



B. Experimental consequences 

If the assumption of incommensurate order is correct 
we would predict a number of interesting results. First, 
there should be four frequencies: the electron pocket 
should be split by the bilayer coupling, although this 
must be small because otherwise it would have been al- 
ready resolved. This sets the bilayer splitting in the 
range of 10 — 20 meV at the level of 10% doping. In 
the overdoped regime the experimentally observed bi- 
layer splitting of 88 meV is already known to be sub- 
stantially smaller than the estimate from the electronic 
structure calculations, which is about 300 meV.'^ Why 
bilayer splitting is so small compared to electronic struc- 
ture calculations is an important theoretical question re- 
quiring attention. To satisfy Luttinger sum rule, we not 
only predict a big hole pocket, but also a smaller hole 
pocket of frequency 250 T. Although there is a very pre- 
liminary indicationilSlit needs to be investigated in further 
measurements and may constitute a confirming evidence 
of either a spiral SDW or an incommensurate DDW. 

IC-DDW could in principle be observed in Y-NMR 
measurements, as the fact that Y atom is situated at 
a high symmetry point is no longer an issue. However, 
the magnetic field at this site arising from the IC-DDW 
currents may be too small to observe this splitting. 

The most dramatic prediction of IC-DDW is the non- 
reflection symmetric band structure in a given spin sec- 
tor. A direct observation of such a band structure will be 
possible in a spin resolved ARPES measurement. Note 
that for spiral SDW up and down spins are mixed. Neu- 
tron scattering intensities from such quasiparticle bands 
may not be sufficient to detect it. 



14 



C. Unresolved issues and future directions 

On the theoretical side a firm grasp on the existence 
of IC-DDW is necessary. This is beyond the scope of the 
present paper whose fundamental concern was to explore 
the phenomenological consequences of IC-DDW. As it 
stands, theory seems to favor the commensurate picture, 
although the predicted hole-pocket in the commensurate 
case is yet to be observed in quantum oscillation measure- 
ments. We expect that many of the experimental puzzles 
will be resolved in future experiments. It is of course de- 
sirable to go beyond Hartree-Fock to formulate a correct 
microscopic theory of IC-DDW. In particular, diagonal 
domain walls as a mechanism for incommensuration have 
not yet been explored. This seems necessary because it 
seems doubtful that the antiphase domain walls that re- 
sult in higher order commensurability, be it from DDW 
or antiphase spin stripe, can explain the present exper- 
iments. Similarly, a rigorous theory of spiral SDW is 
desirable. As to IC-DDW, it is necessary to go beyond 
the approximation of retaining only the single harmonic 
and to formulate a theory that fully conserves currents 
at the vertices. 

The chief experimental puzzle appears to be the con- 
flict of the picture of the Fermi surface emerging from 
the q uantu m oscillations with the ARPES measure- 
ments.'^^EH Although the hole pocket is not a major con- 
cern, because it has been argued that the back side of 
the hole pocket may not be visible because of the DDW 
coherence factors^*^' and similar coherence factors in other 
forms of competing mechanism, the non-existence of the 
electron pockets is more of a concern. From our previous 
calculations it is easy to see that at least two sides of the 
electron pockets should survive the effects of the coher- 
ence factors. It is noteworthy that in the past ARPES 
was able to detect both electron and hole pockets in elec- 
tron doped cuprates.l^ 

Although we do not believe that the observed nega- 
tive Hall coefficient, its magnitude, and its oscillations 
can arise merely from the vortices in the mixed state, be- 
cause it is difficult in such a scenario to obtain multiple 
Fermi pockets tightly constrained by the Luttinger sum 
rule, but it would be interesting to explore further the 
oscillations of the the Hall coefficient, in fact it seems 
almost imperative. 

We have arguecP^ previously that the ortho-II poten- 
tial of YBa2Cu306.5 is not an important factor in the ex- 
planation of the oscillations in contrast to the discussion 
in Ref. |45| . The highly polarizable BaO-layers next to 
the chains should screen the potentials quite effectively, 
to the extent that even disordered chains appear to have 



little effect in the planar physics in many properties. It 
is important to recognize that SdH, and dHvA,E2] mea- 
surements are also available in YBa2Cu408 (~ 14% dop- 
ingEQ), a stoichiometric compound with intrinsic oxygen 
order, where even the negative Hall coefficient is clearly 
observed, in agreement with the observations in ortho-II 
materials. ^^ This degree of universality is unachievable if 
the chain potentials were playing an important role. 

A number of competing mechanismgiSI^ZM] ^q explain 
these unusual experiments have been proposed and have 
been discussed in a previous work by one of us.^^ Thus, 
there is no need to duplicate the discussion here, but the 
salient questions that must be explained are worth re- 
peating here. Are the oscillations, the negative sign, and 
the magnitude of the Hall coefficient explainable with- 
out invoking a two-band scenario of a hole pocket and 
an electron pocket? Can the experimentally observed 
frequencies and the Luttinger sum rule be correctly re- 
produced? If the spiral spin density wave order or the 
antiphase IC-DDW order is the explanation what would 
be the defining experimental predictions that could be 
tested in further experiments? If the quasipartcles of 
a Fermi liquid are truly responsible for the oscillation 
measurements, what is the role of Mott physics in the 
cuprates? 

The lack of any definitive measurements vindicating 
the static order parameters necessary to explain Fermi 
surface reconstruction (recall that dHvA is an equilib- 
rium effect) is worth exploring. A magnetic field even 
as large as 60 T is not sufficient to provide enough per- 
turbation to energetically nucleate static order in YBCO 
with an order parameter large enough to reproduce the 
measurements, especially in YBa2Cu408 with about 14% 
doping. We hope that as other cuprates are explored and 
the theoretical tools are sharpened, many of these ques- 
tions will be answered and will provide a resolution of 
the enigma of high-T^. 
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